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> ■ Abstract 
O 

\& ■ A review of neutron form factor measurements is given. We focus on 

recent measurements of the neutron magnetic form factor G^, and dis- 
■^j- . cuss in detail our measurements of this quantity at momentum transfers 

o : 

O. Q 2 = 0.1 - 0.6 (GeV/c) 2 . 



I. NUCLEON FORM FACTORS 



In the one photon exchange approximation, the elastic electron - nucleon differential 
cross section is described by the so-called Rosenbluth formula@. The expression contains 
two electromagnetic form factors which depend on the four momentum transfer squared 
q 2 = —Q 2 in the reaction, only. These form factors, usually called Ge(Q 2 ) (electric) and 
Gm{Q 2 ) (magnetic), can be thought of as the Fourier representation of the electric charge 
and current distributions in the nucleon. This relation, however, does not hold at realistic 
Q 2 , where the recoil of the nucleon and relativistic kinematics are non- negligible. In the Breit 
(or brick wall) frame, which is the kinematical frame in which no energy is transferred, Ge 
and Gm are proportional to the nucleon charge and current matrix elements, respectively. 

In the late 1950s, Hofstadter and collaborators initiated a series of experiments in which 
the form factors of the proton and the neutron were measured!!. The proton form factors 
were obtained from the measurement of the elastic electron-proton (e — p) cross sections at 
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two different kinematics, keeping Q 2 fixed. Since no free neutron targets are available, the 
neutron form factors were obtained by measuring the quasi-elastic electron-deuteron (e — d) 
cross sections at the same kinematics, and subsequently subtracting the contribution from 
the proton using the e — p data. 

Most measurements of the proton and neutron form factors, even recent onesiHl, have 
used this very same technique. By now, the proton form factors are conventionally considered 
to be known accurately: The proton magnetic form factor G P M is known with an accuracy 
of 2 % up to at least Q 2 = 1 (GeV/c) 2 . The electric form factor G P E is known to about 5 % 
at best, due to the fact that its numerical value is smaller than G P M in the accessible range 
of momentum transfers. In the most elaborate recent work on neutron form factors from 
inclusive measurements, Lung et a/.i have shown that the neutron magnetic form factor G^ d 
can be extracted with an accuracy of 10 - 15 % using this method. The neutron electric 
form factor G\ is essentially undetermined, and recent studies!! suggest that polarization 
variables might be more useful for an unambiguous determination of this small quantity. 

Several measurements have attempted an accurate determination of G\ by employing 
coincident neutron detection, initiated in 1962 by the measurement of Stein et a& Neutron 
detection involves a dedicated neutron detector, usually built from scintillating material, and 
a careful calibration of its detection efficiency Because of these experimentally non-trivial 
considerations, the early measurements^!^ did not yield both reliable and accurate data, and 
consequently large systematic uncertainties had to be introduced. Specific problems were 
the position dependence of the neutron detection efficiency, which was essentially unknown!, 
and the large proton contribution which forced the experimenters to use heavy shieldingjUill 
or to measure only at high momentum transfers^. 

With the advent of high duty cycle electron machines (ELSA (Bonn), MAMI (Mainz), 
Jefferson Laboratory (Newport News)) accurate determinations of the neutron form factors 
are considered realistic, which has lead to a renewed interest in these quantities. 
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II. DETERMINATION OF THE NEUTRON MAGNETIC FORM FACTOR 



We have determined the magnetic form factor of the neutron, G^ f , at four different 
momentum transfers between Q 2 = 0.1 and 0.6 (GeV/c) 2 . The goal of the measurement was 
to minimize the experimental and theoretical uncertainties and to deliver the first accurate 
determination of the behavior of G\ f over an appreciable range of Q 2 . 

The ratio of neutron and proton yields at quasifree kinematics was measured for the 
reactions 2 H(e,e'n) and 2 H(e,e'p) at momentum transfers Q 2 =0.125, 0.255, 0.417 and 
0.605 (GeV/c) 2 (Kin. I, II, III and IV), detecting the neutron and the proton simultane- 
ously in the same scintillator array, the nucleon detector. The neutron detection efficiency 
was measured in situ with the 1 H(7,7r + )n reaction. From this, the ratio R of 2 H(e,e'n) and 
2 H(e,e'p) cross sections was determined and used to extract the neutron magnetic form fac- 
tor Gyi in a model insensitive approach, resulting in an inaccuracy between 2.1 and 3.3 % 
in G n M . 

The method employed, first suggested by Durandll!, was previously applied in only a 
few experimentsiH. In this way the luminosity, the electron detection efficiency, and the 
electron solid angle cancel out, while the nucleon acceptances and the choice of the deuteron 
wave function cancel out to first order. The lay-out of the experiment minimized corrections 
due to proton losses, and the restriction to quasifree kinematics minimized nuclear effects 
and their corresponding uncertainties. 

The ratio R of the cross sections of the 2 H(e,e'n) and 2 H(e,e'p) reactions was determined 
after taking into account proton losses due to nuclear reactions, multiple scattering and 
edge effects, and an in situ calibration of the neutron detection efficiency. G 1 ^ was extracted 
taking into account nuclear effects, using the known e — p cross sections and available infor- 
mation on G"%. G% contributes less than 2 % to the e — n cross section and is not a dominant 
source of uncertainty. 



3 



A. Experimental details 

1. The detector setup 

The electron beam (between 900 and 1600 MeV, 20 to 60 nA, 20 to 50 % duty cycle, 
depending on kinematics), delivered by the ELS A accelerator at the Physics Institute at 
Bonn, impinged axially on a cylindrical target (length 10 cm) made from 125 //m thick kapton 
and filled with liquid Hydrogen or Deuterium. Electrons (and pions from the 1 H(7,7r + )n 
reaction) were detected in the ELAN magnetic spectrometer with four scintillators and 
MWPCs and one Cerenkov counter. 

Protons and neutrons were detected in the nucleon detector using the time of flight 
method. Since protons are about hundred times more abundant than neutrons, a reliable 
particle identification is important. We used three thin veto scintillators which allowed a 
misidentification rate of far less than 10~ 4 , and the verification of the losses due to pile -upi. 

For Kinematics I, II and III, the nucleon detector, which was used previously, consisted 
of a total of five scintillators (NE102a) of dimensions 25 x 25 cm 2 with a thickness of 
2 mm for the AE detectors (1 mm for Kin. I) and 50 mm for the E counters. The first E 
counter (Ef ront ) detected simultaneously protons and neutrons, while shielding behind Ef ront 
prevented protons from the quasielastic 2 H(e,e'p) reaction to reach Eback- The neutrons 
detected in Eback were used to determine the losses in the neutron yield in Ef ront caused 
by the software gates on the veto detectors. For Kin. IV, the detectors of dimensions 
100 x 18 cm 2 and thickness 1 cm and 18 cm, respectively, had a two-sided readout. At 
this large proton energy the use of lead to shield Eback becomes impracticable. We therefore 
installed two veto detectors between Ef ront and Eback in order to determine the neutron 
losses. For a more detailed description, see RefH Each measurement of R was bracketed 
by two calibration runs for the neutron detection efficiency To switch from the calibration 
to the measurement of the ratio, only the target and the spectrometer settings needed to be 
changed, which reduced the risk of systematic errors. 
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B. The analysis 



Neutrons detected in the nucleon detector were defined by a signal in Ef ront or in E back 
at the correct time with respect to the spectrometer and in excess of a certain software 
adjustable threshold, together with a neutron condition on the veto counters. Different veto 
conditions on the ADC and TDC information were applied and used for stringent cross 
checks. Thanks to the duty cycle of the ELSA stretcher ring (20 to 50 %), a signal to noise 
ratio of 200:1 for neutrons was achieved (Fig. [I]). 

The stability of the neutron detection efficiency of Ef ront , rj n , was monitored by means 
of the ADC signal for protons, and by the number of neutrons coming from the 2 H(e,e'n) 
reaction scaled with the number of electrons in the spectrometer. At Kin. Ill and IV, the 
2 H(7,p)n reaction served as an additional monitor for the stability of rj n . This reaction can 
not be used to determine the absolute value of r) n since pion production cannot be excluded. 

The reaction 1 H(7,7r + )n was used to obtain an in situ, absolute determination of the 
detector efficiency for neutrons. Under the assumption that only real photons are used, 
the kinematics of this reaction can be arranged in such a way as to unambiguously tag 
neutrons impinging on the nucleon detector centered at the same energy as those from the 
2 H(e,e'n) reaction. Eighty percent of the pion flux in the spectrometer was determined to 
actually originate from virtual photons, i.e. from the 1 H(e,7r + )e'X reaction. Monte Carlo 
simulations^!, however, show that more than 99 % of the pions which enter the spectrometer 
(and satisfy the cuts used in the analysis, as described in Ref.0) originate from electrons 
scattering to extremely small angles (<C 5°). Because of the stringent requirements put 
on the pion kinematics, to resulting data set can be treated as originating from the real 
photon induced 1 H(7,7r + )n reaction. This consideration is unfortunately not supported by 
simulations performed by the Basel group, which have recently published competing G 1 ^ 
results^. The discrepancy has to be resolved in the near future. Curiously, for neutrons 
with an average energy of 61 MeV, the measured detection efficiency in the center of the 
nucleon detector as a function of threshold agreed with the one established earlier0'0. 
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In order to obtain R, the measured yields must be corrected for: i) the net proton losses 
due to nuclear reactions and multiple scattering in the material between the reaction vertex 
and the detector; ii) the dependence of the neutron detection efficiency on the distribution of 
the neutrons in space and energy over the detector surface. Details about the small correc- 
tions for Hydrogen contamination of the Deuterium, and for the contributions of the target 
end caps to the different reactions can be found in Ref.0. In addition, for an extraction 
of G^j one must evaluate nuclear effects, such as final state interactions (FSI), meson ex- 
change currents (MEC) and isobar currents (IC), which alter the proton and neutron yields 
expected from free proton respectively neutron targets. 

The 1 H(e,e')p reaction was used for the determination of the proton detection efficiency, 
including the above mentioned losses. The experimentally found losses were consistent with 
numerical checks done with the GEANT packag which was extended to include total 
proton cross sections at low energies^. Protons from the 1 H(e,e')p and the 2 H(e,e'p)n 
reaction also were used to calibrate and monitor the light response of the scintillators in 
view of the threshold dependence of the neutron detection efficiency. 

The reactions 2 H(e,e'n) and 1 H(7,7r + )n lead to different energy and position distributions 
of neutrons in Ef ront . While the neutron distribution for the 1 H(7,7r + )n reaction is known 
from experiment, the 2 H(e,e'n) distribution was obtained using the ENIGMA Monte Carlo 
CO de§. The second order effect of the different distributions (including edge effects) on the 
ratio R was simulated with the program KSUVAX0. The response of the detector used for 
Kin. I, II and III had been calibrated carefully previously^ using a tagged neutron beam at 
PSI. The positional dependence of the response of the detector used for Kin. IV was studied 
by means of its double read out system. 

C. Nuclear effects 

Nuclear effects cause R to differ from the expected ratio Rf ree for free nucleons: Rf ree = 
R*(l — 5R). The corrections SR were calculated in two different models, of which the firstH 
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includes FSI, MEC and IC, and shows that FSI dominate the other corrections, whereas the 
second oneH includes only FSI but is fully relativistic. It appears that the result for G\ f 
obtained from R is affected to less than 0.4 % by relativistic effects. We use the difference 
between the two calculations, including FSI only, as an indication of the inherent model 
uncertainties, and their average as the most probable correct value. To this average, we add 
then the effects of MEC and IC which are taken as the difference between the full calculation, 
and the calculation including FSI only, both within the same mo deli. We conservatively 
estimate the uncertainty due to these effects to be half their size. The values of SR range 
from 2.2 to 8.5 %, depending on the kinematics. 



D. Results 

For the extraction of G^, the e — p cross sections are needed. The measured cross 
section dataH have been averaged in small intervals around the central Q 2 values, tak- 
ing into account the Q 2 dependence obtained from the dipole form factors for the proton 
(G p E =G D =(l+Q 2 /a 2 )- 2 , G P M =2.793*G P E , and a 2 =0.710 (GeV/c) 2 ). The statistical errors of 
the used data points determine the statistical error in the average. In order to obtain the sys- 
tematical error, the average has been recalculated after each e—p data set had been changed 
individually by its systematical errors, and adding the resulting changes quadratically. G% 
was chosen to be 0.037 ± 0.0170. 

The results for G\ are shown in Fig. ^|. We confirm our pilot measurement at 
Q 2 =0.125 (GeV/c) 2 , performed at the MEA accelerator (NIKHEF, Amsterdam) with the 
same nucleon detector as used in Kin. I, II and III. Note that the published results differ by 
more than one standard deviation, only because of the difference in the value for the proton 
cross section used. This difference, 2 %, comes from the different Q 2 range of the e — p data 
used in the fit and reflects our knowledge of the e—p cross section. Our value for G 7 ^ at this 
momentum transfer is significantly higher than our previous result for Q 2 =0.093 (GeV/c) 2 
obtained at MEA. Fig. |2] suggests a slope in G^/Gd around these momentum transfers. 
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Recently Markowitz et al.EJJ determined G 1 ^ from an absolute measurement of the 
2 H(e,e'n) cross section. At 0.109 (GeV/c) 2 , their result lies 13.6 % above the dipole fit, 
whereas our result at 0.125 (GeV/c) 2 is only 2.5 % above. Relativistic effects, which have 
not been applied in Ref.0, do not appear to remove this discrepancy. Our results agree with 
the other results of Markowitz, and with the result obtained by Gao@ using polarisation 
techniques. Fig. |2| shows that most of the recent theoretical predictions of G n M are not in 
agreement with our data. We only show some of the (rare) theoretical predictions which 
are not fits to data. The only model which reproduces the trend indicated by the data 
is the minimal model of Meissneril. Since Meissner's complete model does not reproduce 
the suggested trend, we observe that none of the parameter-free models presently available 
describes the trend in G\ sl as observed in our data. The data show for the first time a 
significant violation of the scaling rule G p M /fi p = G^/ '//„. 



III. CONCLUSIONS AND OUTLOOK 

We have combined the ratio method proposed long ago with the possibilities offered only 
recently by the large duty cycle electron beams to measure with substantially decreased 
uncertainties. The coincident detection of the electron and the knocked out nucleon reduces 
significantly uncertainties due to nuclear effects, whereas the large duty cycle allows to detect 
simultaneously protons and neutrons. 

Recently, a new set of data has become available from the MAMI accelerator, MainzEH 
covering a similar range in momentum transfer. Although these data confirm the violation 
of the above mentioned scaling rule, they differ significantly from our results. This fact is 
under investigation, and a possible explanation has been given in section IIB. 

An experiment has been approved at the Bates accelerator, in which G\ will be measured 
at small Q 2 , i.e. < 0.125 (GeV/c) 2 , using the method as employed by Ref.0. Although 
nuclear effects are large for the proposed measurement, the result will be a continuous set 
of data points up to the region where the existence of a sudden fluctuation in G 1 ^ has been 
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claimecfi 

At TJNAF, an experiment has been approved in which G 1 ^ will be measured using the 
ratio method at Q 2 = 0.3 — 5.1 (GeV/c) 2 . One of the main experimental problems in this 
experiment will be the matching of the proton and neutron phase spaces, due to the different 
techniques employed to detect each of the nucleons. 
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FIG. 1. Time spectra of protons (line) and neutrons (shaded area) using a subset of the data 
from the measurement at Q 2 =0.255 (GeV/c) 2 . 
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FIG. 2. Data from recent G"%j measurements up to Q 2 =0.8 (GeV/c) 2 , scaled to the dipole fit: 
black circles [17], open squares [11], open circle [25]. The black squares are from this work. Several 
model predictions are shown. VMD: The minimal (dotted line) and complete (long dash) model 
of [28]. Constituent Quark Model: non-relativistic calculation of [29] (dot-long dash), relativistic 
calculation of [30] (dot-short dash). 
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